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We calculate the photoproduction of heavy quarks in proton- proton collisions at RHIC, Tevatron 
and LHC energies, where the photon reaches energies larger than those ones accessible at DESY- 
HERA. The integrated cross section and the rapidity distributions for open charm and bottom 
production are computed employing sound high energy QCD formalisms. For the linear pQCD 
approaches we consider both the usual collinear factorization and the fcx-factorization formalisms, 
■ whereas for the nonlinear QCD (saturation) calculations one considers the Golec-Biernat-Wvisthoff 

and the Iancu-Itakura-Munier parameterizations for the dipole cross section within the color dipole 
f^*) , picture. 
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INTRODUCTION 



Heavy quark production in hard collisions of hadrons, leptons, and photons has been considered as a clean test of 
perturbative QCD (For a recent review see, e.g., Ref 0). This process provides not only many tests of perturbative 
QCD, but also some of the most important backgrounds to new physics processes, which have motivated an extense 
phenomenology at DESY-HERA, Tevatron and LHC. These studies are mainly motivated by the strong dependence of 
the cross section on the behavior of the gluon distribution, which determines the QCD dynamics at high energies. In 
particular, the heavy quark photoproduction on nucleon and nuclei targets has been studied in detail in Refs. |jj>|^,Q, 
considering the several available scenarios for the QCD dynamics at high energies. The results of those analysis show 
that future electron-proton (nucleus) colliders at HERA and RHIC PJlIj probably could determine whether parton 
distributions saturate and a stringent constraint to the behavior of the gluon distribution in the full kinematical range 
Q_i| could be posed. Along these lines, in Ref. we have analyzed the possibility of using ultraperipheral heavy ion 
JL . collisions (UPC's) as a photonuclear collider and studied the heavy quark production assuming distinct formalisms 
for the QCD evolution. 

Recently, Klein and Nystrand |8( have analyzed the quarkonium photoproduction in proton-proton collisions, con- 
sidering the energetic protons as search of large electromagnetic fields. In particular, these authors have used the 
. t-H ! photon spectrum from Ref. and a photon-proton cross section for quarkonium production obtained by fitting the 
HI and ZEUS data. The main conclusion of that study is that the cross sections are large enough for this channel 
to be observed experimentally and that this reaction can be used to study the gluon distribution in protons at small 
values of the Bjorken x variable. This achievement motivates the analysis of other processes that are sensitive to the 
gluon distribution and have larger cross sections. 

In this paper, we study the photoproduction of heavy quarks in proton-proton collisions considering distinct theo- 
retical scenarios, which have been analyzed in detail in our previous papers 0,0,0- Our main motivation comes from 
the fact that in this process the photon reaches energies higher than those currently accessible at DESY - HERA. 
Similar motivation is present in Ref. 0, where vector meson photoproduction is investigated. Here, we estimate, 
for the first time, the total cross section and the rapidity dependence of the photoproduction of heavy quarks in 
proton-proton collisions, considering the fc^-factorization approach taking into account distinct unintegrated gluon 
distributions. Moreover, based on the color dipole picture we compute the total cross section for charm and 
bottom photoproduction within the Iancu-Itakura-Munier (IIM) model for the dipole cross section [l^ . These results 
will be used as input for our further calculations in proton-proton collisions. For comparison, we also present the 
predictions from the collinear factorization approach. 

In relativistic heavy ion colliders, the heavy nuclei give rise to strong electromagnetic fields, which can interact with 
each other. In a similar way, these processes also occur when considering energetic protons in pp(p) colliders. Namely, 
quasi-real photons scatters off protons at very high energies in the current hadron colliders. In particular, the heavy 
quark photoproduction cross section in a proton-proton collision is given by, 

a(p + p^p + X + QQ) = 2^ ^Ma^ Q Q X (w^ p = 2toy/S^) du , (1) 
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where u> is the photon energy in the center-of-mass frame (c.m.s.), W lp is the c.m.s. photon-proton energy and \/Snn 
denotes the proton-proton c.m.s. energy. The photon spectrum is given by 0, 

A ^ 11 3 3 1 \ 

with the notation 0=1 + [ (0.71 GeV 2 VQ^J and Q^ in = w 2 /^ (1 - 2u/y/S NN )] w {uj/il) 2 , where 7 L is the 
Lorentz factor. The expression above is derived considering the Weizsacker- Williams method of virtual photons and 
using an elastic proton form factor (For more detail see Refs. 0,0)- It is important to emphasize that the expression 
PJl is based on a heuristic approximation, which leads to an overestimation of the cross section at high energies ( 

11% at i/s = 1.3 TeV) in comparison with the more rigorous derivation of the photon spectrum for elastic scattering 
on protons derived in Ref. .10]. For a more detailed comparison among the different photon spectra see Ref. [Tlj . 

Other process of interest for heavy quark production is the coherent interaction between two photons. As verified 
in Ref. |12j] for ultraperipheral heavy ion collisions, the QCD dynamics implies an enhancement of the cross section 
in comparison with Quark - Parton Model calculations |L|. However, its cross section is a factor about 800 smaller 
than for the one - photon process. Furthermore, differently from two-photon interactions, which leave both photon 
intact, photoproduction should only dissociate the proton target. Then one way to select photoproduction events is 
to eliminate events where both protons breakup and requiring a single rapidity gap. 

In what follows, we briefly present sound models for the heavy quark photoproduction at the photon level, a tot (jp — > 
QQX), which enter as input in the calculations for proton- proton collisions in Eq. (JTJ. Our results are compared 
with the DESY-HERA data 0] and the prediction from the IIM model for this process is presented for the first 
time. In the last section, the numerical results for the rapidity y of the produced states and their total cross sections 
are shown. A comparison on the order of magnitude of the cross sections for the distinct approaches is performed. 
Moreover, we present our main conclusions. 
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II. A COMPARISON AMONG HIGH ENERGY APPROACHES 

The photon-proton cross section can be calculated considering different theoretical scenarios Usually, one 

calculates it assuming the validity of the collinear factorization, where the cross sections involving incoming hadrons 
are given, at all orders, by the convolution of intrinsically non-perturbative, but universal, quantities - the parton 
densities, with perturbatively calculable hard matrix elements, which are process dependent. In this approach, all 
partons involved are assumed to be on mass shell, carrying only longitudinal momenta, and their transverse momenta 
are neglected in the QCD matrix elements. The heavy quark cross section is given in terms of the convolution between 
the elementary cross section for the subprocess 75 — > QQ and the probability of finding a gluon inside the proton, 
namely the gluon distribution. The photoproduction cross section at leading order is given by [18|. 



atot{lP^ QQX) = 47re Q / 2 W a em a s (^ F ) g p (x, fip) 



QQ 



1 + ^H'°(t ± ^)- (1 + /,) ^ 



(3) 



where Mqq is the invariant mass of the heavy quark pair, with x = M^^/W^p and g P (x,fj,p) is the gluon density 

inside the proton at the factorization scale [Pp. In addition, mg is the heavy quark mass, eg is its electric charge and 
(3 = 4 mg/ 'Mqq- For the present purpose we will use [P F = imq, with m c = 1.5 GeV and mj, = 4.5 GeV. In our further 
calculations on the collinear approach one takes the gluon distribution given by the GRV98(LO) parameterization 
[H^|. It should be noticed that different choices for the factorization scale and quark mass produce distinct overall 
normalization to the total cross section at photon-nuclcon interactions and that NLO corrections can be absorbed in 
these redefinitions of [P F and mg. In the next section we discuss in more detail the dependence of our results in the 
parton distributions used as input in our calculations. 

On the other hand, in the large energy (small- 2; ) limit, the effects of the finite transverse momenta of the incom- 
ing partons become important, and the factorization must be generalized, implying that the cross sections are now 
/c^-factorized into an off-shell partonic cross section and a fc^-unintegrated parton density function Fix, k±), charac- 
terizing the fc^-factorization approach [2(|E],|22. The function T is obtained as a solution of the evolution equations 
associated to the dynamics that governs the QCD at high energies (see ^3 f° r a review). The latter can recover 
the usual parton distributions in the double leading logarithmic limit (DLL) by its integration over the transverse 
momentum of the fcj_ exchanged gluon. The gluon longitudinal momentum fraction is related to the c.m.s. energy, 
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W 1P , in the heavy quark photoproduction case as x = A.rn,Q/W 2 p . This assumption is a very good approximation, 
thought the scaling variable x in fact depends on the kinematic variables of the incoming particles, namely on z and 
parton momenta (for details on these issues see 2] and references therein). The cross section for the heavy-quark 
photoproduction process is given by the convolution of the unintegrated gluon function with the off-shell matrix ele- 
ments 0, UJl HE Hj| • Considering only the direct component of the photon we have that the total cross section reads 
as 0, 

e o /" d 2 k± 

7T J fej_ 

where D\ = p 2 ± + rrig and D2 = (fej_ — P\±) 2 + ttVq. The transverse momenta of the heavy quark (antiquark) are 
denoted by p 11 _ and p 2 j_ = (^JL — Pxx)i respectively. The heavy quark longitudinal momentum fraction is labeled by 
z. The hard scale /i in general is taken to be equal to the gluon virtuality, in close connection with the BLM scheme 
|25j . Therefore, for our purpose we will use the prescription /x 2 = k\ + Trig. As some evolution equations for the 
unintegrated gluon distribution, for instance the CCFM equation j2||, enable also its dependence on the scale fj,, we 
have explicitly written down it in Eq. J3J for sake of generality. 

In performing a numerical analysis within the k± -factorization approach (hereafter labeled SEMIHARD), we use 
two distinct parameterizations for the unintegrated gluon distribution (For details see Ref. p|). First, one considers 
the derivative of the collinear gluon distribution, quite successful in the proton case Q and investigated in the nuclear 
case in Ref. 0. It simply reads as 

•'DLL [X, fcjj = d lnfc 2 > ( 5 ) 

where xg(x, Q 2 ) is the gluon distribution, which was taken from the GRV98(LO) parton distribution [l9{. A shortcom- 
ing of the function above is that it eventually produces negative values for the unintegrated gluon distribution at large 
x. This has consequences in the description of the behavior near threshold, being more important for bottom pro- 
duction than for charm. This is due to the scaling variable to be proportional to the quark mass, x = Araq/W 2 p , and 
in general the calculation slightly underestimates the cross section value at low energies. A more accurate calculation 
consists on multiplying the collinear gluon function above by the Sudakov-like form factor T g (k, /j) [27j . 

As the heavy quark production is characterized by a scale of order of the quark mass, the charm production, in 
particular, is considered a primary candidate for investigating the kinematic region of QCD on the boundary between 
perturbative and high density QCD (For recent reviews see, e.g., Refs. |2S,|2g). Recently, we have demonstrated 
that the inclusive charm total cross section exhibits the property of geometric scaling [3(j , which is one of the main 
characteristics of the high density approaches. This fact motivates to estimate these features in the photoproduction 
case. In order to do so, we consider the Golec-Biernat-Wiisthoff parameterization [3l| (hereafter SAT-MOD), which 
gives the following unintegrated gluon distribution, 



where Q 2 at ( x ) is the saturation scale, which defines the onset of the nonlinear QCD (saturation) effects. The param- 
eters were obtained from a fit of the inclusive structure function F2 and total photoproduction cross section data at 
DESY-HERA data. The fit using the expression of the color dipole picture, as explained below, and the saturation 
ansatz for the dipole cross section produced Co = 29.12 mb, A = 0.277 and Xq = 0.41 ■ 10~ 4 . This fit procedure also 
includes the charm contribution. In addition, we have included the large- x threshold factor (1 — x) 7 (for the heavy 
quark case) in Eq. © to take into account the correct energy dependence at low energies near production threshold. 

The saturation model is based on the color dipole picture of the photon-proton interaction |14j . In the proton rest 
frame, the DIS process can be seen as a succession in time of three factorisable subprocesses: i) the photon fluctuates 
in a quark- antiquark pair with transverse separation r ~ 1/Q long after the interaction, ii) this color dipole interacts 
with the proton target, iii^ the quark pair annihilates in a virtual photon. The interaction j*p is further factorized in 
the simple formulation \ l 



.7 p 

T L,T 



(x, Q 2 ) = J2 / dz d2r I*£,t ( 2 . r > Q 2 )\ 2 °d ip (x, r) , (7) 
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where z is the longitudinal momentum fraction of the quark. The photon wavefunctions L T are determined from 
light cone perturbation theory and read as 

l^l 2 = ^4 {[^ + (1 - *)V*!(*r) + m) K 2 {er)} 

\H\ 2 =^pe}{Q 2 zHl-zfKUer)}, (8) 

where the auxiliary variable s 2 = z(l — z) Q 2 + m 2 depends on the quark mass, m/. The -Ko.i & r e the McDonald 
functions and the summation is performed over the quark flavors. The dipole-hadron cross section a^n contains all 
information about the target and the strong interaction physics. In general, the saturation models [ljj, I2J |3^, |3^| 
interpolate between the small and large dipole configurations, providing color transparency behavior, a dip ~ r 2 , as 
r Q sa t and constant behavior at large dipole separations r < Qsat- Along these lines, the phenomenological 
saturation model resembles the mains features of the Glauber-Mueller resummation. Namely, the dipole cross section 
takes the eikonal-like form, 



crdipix, r) = a Q 



1 — exp 



(9) 



Its phenomenological application has been successful in a wide class of processes with a photon probe. An important 
aspect of the saturation models is that they resume a class of higher twist contributions which should be non-negligible 
towards the low Q 2 region 0, • Moreover, it is important to emphasize that the dipole approach was extended for 
heavy quark production in proton-proton (nucleus) collisions in Refs. l36|. [37ll. and its equivalence with the collinear 
approach as well as the comparison with pp data was presented in Ref . 38] . 

Although the saturation model describes reasonably well the HERA data, its functional form is only an approx- 
imation of the theoretical non-linear QCD approaches. On the other hand, an analytical expression for the dipole 
cross section can be obtained within the BFKL formalism. Currently, intense theoretical studies has been performed 
towards an understanding of the BFKL approach in the border of the saturation region [39l l40j . In particular, the 
dipole cross section has been calculated in both LO and NLO BFKL approach in the geometric scaling region |41| . It 
reads as, 



o-di P (x,r) =<7 Q [r 2 Q 2 at (x)] sat exp 



In 2 (r 2 QL 

2(3a s Y 



(10) 



where ctq = 2nR 2 (R p is the proton radius) is the overall normalization and the power 7 sat is the (BFKL) saddle point 
in the vicinity of the saturation line Q 2 = Q 2 at ( x ) (the anomalous dimension is defined as 7 = 1 — 7 sa t). As usual in 
the BFKL formalism, a s = N c a s /ir, (3 ~ 28C(3) and Y = ln(l/a;). The quadratic diffusion factor in the exponential 
gives rise to the scaling violations. 

The dipole cross section in Eq. (|1U|) does not include an extrapolation from the geometric scaling region to the 
saturation region. This has been recently implemented in Ref. 0, where the dipole amplitude N(x, r) = adi P /2nR 2 
was constructed to smoothly interpol between the limiting behaviors analytically under control: the solution of the 
BFKL equation for small dipole sizes, r <C l/Q s &t{x), and the Levin- Tuchin law 01 f° r larger ones, r ^> 1/Q sa t(^)- 
A fit to the structure function F%{x,Q ) was performed in the kinematical range of interest, showing that it is not 
very sensitive to the details of the interpolation. The dipole cross section was parameterized as follows, 



^p° (*,r) = <?o (2 J 7 , for rQ sat (x) < 2 , (n) 

[ 1 - exp [-a In 2 (brQ sat )] , for rQ sat (a;) > 2 , 

where the expression for rQ s& t{x) > 2 (saturation region) has the correct functional form, as obtained either by 
solving the Balitsky-Kovchegov (BK) equation or from the theory of the Color Glass Condensate (CGC) |44j |. 
Hereafter, we label the model above by CGC. The coefficients a and b are determined from the continuity conditions 
of the dipole cross section at rQ sa ,t(x) = 2. The coefficients 7 sat = 0.63 and k = 9.9 are fixed from their LO BFKL 
values. In our further calculations it will be used the parameters R p = 0.641 fm, A = 0.253, xq = 0.267 x 10 -4 and 
Afo = 0.7, which give the best fit result. We have included also a large-x factor as for the saturation model. In the 
color dipole picture, the heavy quark photoproduction cross section using the CGC model reads as, 

(Ttot (TIP -» QQX) = fdzd 2 r\^(z, r,Q 2 = 0)\ 2 a^ c (x,r), (12) 



where the longitudinal piece does not contribute as I^&lI 2 cx Q 2 . The transverse contribution is computed using Eq. 
© and introducing the appropriated mass and charge of the charm or bottom quark. 
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The dipole approach is very useful in providing semi-analytic solutions for the cross section. In the heavy quark 
case, it can be shown that the scattering process is dominated by small size dipoles with mean value r ~ 1/wq. 
In this scenario, the photon wavefunction selects the color transparency region as the main contribution to the 
cross section. Concerning the longitudinal fraction z of the dipoles, at the limit Q 2 -C 4nig only the symmetric 
configurations <z>« 1/2 contribute. In this case, which it is obeyed in the photoproduction regime, we can write in 
a semi-quantitative way, 



where a is the small-r limit of the dipole cross section. For heavy quark production we have that the saturation model 
predicts a — uq Q 2 at r 2 /4, whereas for the CGC model it reads as a = <x Ao {Q 2 a t r 2 /4) 7aat modulo the diffusion term 
in Eq. (JTJJ. Therefore, by introducing these values in Eq. (fT3f one produces the following analytical results for both 
models, 



which implies in general a slightly smoother energy growth for CGC than for SAT-MOD since 7 sa t < 1. It should 
be noticed that this behavior is obtained only in the case of sufficiently small x (large rapidity Y) , where we could 
neglect the diffusion term for CGC and take only the anomalous dimension at the saturation vicinity. On the 
other hand, for not so small x the diffusion factor can not be disregarded and the effective anomalous dimension 
7eg (x, r) = 7sat + [ m (2/rQ sat )/K A Y ] in Eq. 1)1 1H should be taken. This situation can occur in the bottom production 
at intermediate energies due to the large bottom mass, which implies the CGC having practically the same slope on 
energy as the saturation model. This can be understood in realizing that the effective anomalous dimension for 

Q A rr-\ ~\J[ f~\T~\ 

SAT-MOD is 7° ff ~ =1 and at sufficiently large x they are very close of each other. 

Having presented the main sound approaches for heavy quark photoproduction at the photon level, let us compare 
their numerical results with the experimental DESY-HERA data [lfj. They are shown in Fig. ^ where the following 
notation is considered: the results from the saturation model (SAT-MOD) are denoted by the solid lines and the usual 
collincar factorization calculation is labeled by the dot-dashed curves. The fcj_-factorization formalism (SEMIHARD), 
using the unintegrated gluon function in Eq. (J5J) produces the long-dashed curves, where the Iancu-Itakura-Municr 
parameterization for the dipole cross section (CGC) gives the dashed curves. In all calculations we have use the same 
quark masses m c = 1.5 GeV and m& = 4.5 GeV. We quote Refs. @>E1 for detailed investigations on the dependence 
of the cross section in choosing different collinear gluon parameterizations and different factorization scale for both 
the collinear and fcj_-factorization approaches. 

Lets compare the distinct behaviors on energy and overall normalizations. Both the collinear approach and the 
semihard formalism give a consistent data description at high energies. At the region near the threshold, the semihard 
approach slightly underestimate the cross section, being more pronounced in the bottom case. This is understood 
in terms of the considered unintegrated gluon distribution J-dll, as discussed before. The saturation model under- 
estimates the cross section at high energies by a factor about 2, producing a reasonable description of the region 
near threshold. The CGC model follows similar trends, but the high energy values are closer to the experimental 
measurements than the saturation model. The slope on energy for SAT-MOD and CGC can be easily understood 
from the semi-quantitative results in Eq. I|14fl . For the bottom case, we are not in sufficiently small x and SAT-MOD 
and CGC give similar energy slope in agreement with our discussion above. Unfortunately, the current precision and 
statistics of the experimental measurements of the photoproduction cross section are either low to formulate definitive 
conclusions about the robustness of the different approaches presented here. More precise measurements could be 
pose stringent constraints on the energy dependence and overall normalization. Finally, it should be noticed that the 
present calculations concern only the direct photon contribution to the cross section, whereas the resolved component 
has been neglected. In some extent the results from the saturation models presented here let some room for this 
contribution. Details on its calculation and size of its contribution can be found, for instance, in Ref. 24] . 



In what follows, we will compute the rapidity distribution and total cross sections for the photoproduction of open 
charm and bottom from proton-proton collisions at high energies. The approaches shortly reviewed in the previous 
section serve as input for the numerical calculations using Eq. Q for the energies of the current and future pp and 
pp accelerators. Namely, one considers the shorter pp running at RHIC upon energy of \/Snn — 500 GeV and the 




(13) 




(14) 



III. RESULTS AND CONCLUSIONS 
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Tevatron value ^/Snn = 1-96 TeV for its pp running. For the planned LHC pp running (ATLAS and/or CMS) one 
takes the design energy ^/Snn — 14 TeV. 

The distribution on rapidity y of the produced open heavy quark state can be directly computed from Eq. 
by using its relation with the photon energy u>, i.e. y oc ln(w/mg). A reflection around y — takes into account 
the interchanging between the proton's photon emitter and the proton target. Explicitly, the rapidity distribution is 
written down as, 

da [p + p{p) ^ QQ + X + p(p)] JV (a,) 

dy = u — H ■ (15) 

Lets start by analyzing the predictions from the collinear approach, obtained using distinct gluon pdf's. In Fig. 
12 the comparison among the resulting rapidity distributions for open charm photoproduction at LHC energy is 
presented. It has been considered the GRV (94 and 98) and MRST (2001) parton parameterizations. The results for 
bottom and/or charm at RHIC/Tcvatron energies give somewhat very close curves and they will be not presented here. 
The results presented in Fig. [21 are strongly dependent on the gluon distribution, which motivates the study of this 
process in order to constrain it. Moreover, they can useful in understanding what is the underlying QCD dynamics 
at high energies. We have that the deviations between the predictions from the GRV98(LO) and MRST2001(LO) 
pdf's are approximately of 30 %. In order to demonstrate the strong dependence of the rapidity distribution on the 
gluon density, we also present the result obtained from the obsolete GRV94(LO) pdf, which it has a steeper gluon 
distribution at small x. In our further comparisons with other approaches, we will use the GRV98 pdf. 

The resulting cross sections coming out of the distinct theoretical inputs considered in previous section are depicted 
in Figs, (l.'il')ll at RHIC, Tevatron and LHC energies, respectively. Lets discuss in general lines the present results. 
At RHIC energies, larger values of x oc mq/W lp are probed in the process than at Tevatron and LHC, mostly for 
bottom. This fact explains the very similar result for all approaches in the bottom case in Fig. H (right plot). For the 
charm case, the situation starts to be different. The SEMIHARD result presents a sharper tail at large rapidities as a 
consequence of its behavior on energy near threshold, as discussed before. The overall normalization and its behavior 
at central rapidities follow the same trend coming from the cross section at photon level (see Fig. QJ. Collinear and 
/c^-factorization formalisms often give closer results, whereas the saturation models (SAT-MOD and CGC) produce 
smaller cross sections. 

At Tevatron energies (Fig. 0J, smaller x are being probed in the reaction, and the cross section becomes more 
dependent on the high energy behavior of the cross section at photon level. The overall normalization becomes 
increasingly sizeable for different approaches. For the charm case, the result could already distinguish between the 
saturation approaches and the usual pQCD calculations. At LHC energies (Fig. the separation is more clear once 
smaller x are probed. High values are obtained from the collinear approach due to its steeper energy growth, followed 
by the semihard approach. The latter still considers an unintegrated gluon function dependent on the collinear one 
and therefore has closer behavior on energy. Concerning the saturation approaches, SAT-MOD provides the lower 
limit, whereas CGC gives somewhat larger values. This is once again a consequence of their high energy behavior, as 
verified in Fig. ^ 

Let us now compute the integrated cross section considering the distinct QCD approaches. The results are presented 
in Table for the open charm and bottom pair production at RHIC, Tevatron and LHC, respectively. The collinear 
factorization approach gives the largest rates among the models studied, followed by the semihard formalism, as a 
clear trend from the distribution on rapidity. Concerning the saturation models, the CGC results in a cross section 
of order 20 % larger than SAT-MOD. The values are either large at Tevatron and LHC, going from some units of /xbs 
at RHIC to hundreds of /As at LHC. Therefore, these reactions can have high rates at the LHC kinematical regime. 
As stated before, we don't consider the resolved photon contribution in our calculations. However, it is possible to 
present an estimate of this process considering the results from Refs. [l3L l45| for ultraperipheral heavy ion collisions, 
where this contribution was studied using collinear factorization. One of the main results is that these contributions 
are « 15% and 20% of the total charm and bottom photoproduction cross sections at LHC energy, respectively. On 
the other hand, in Ref. [24[ this contribution was estimated to be of order of 20-30 % of the direct photon cross 
section using fcj_-factorization. It is important to emphasize that the inclusion of resolved photon contribution brings 
the predictions for bottom production closer to the jp data, but a major inconsistency cannot be claimed due to the 
large experimental errors and theoretical uncertainties. We postpone a detailed analysis on these issues for a future 
publication. 

The cross section computed here can be contrasted with the open charm hadroproduction in central pp(p) colli- 
sions. Recently, Raufeisen and Peng 38] have computed it considering the NLO parton model and the color dipole 
formulation. They found that the results are subject to uncertainties coming from different choices for quark mass 
and parameters of the models, but they are able to describe all data (except recent STAR measurement) inside the 
uncertainty band. The photoproduction cross section is of order 0.1 % from the corresponding hadroproduction. 
However, as we will see below, the production rates are either high and the experimental signal is significantly clear. 
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At RHIC, where the luminosity is assumed to be £rhic = 10 31 cm _2 s -1 , the open charm rate ranges on 3.7— 7.8x 10 6 
events by year. It should be noticed that RHIC uses most part of the running in the heavy-ion mode, so we used 10 6 
s in the last estimation. For the bottom case, the rates are 3.7 — 4.3 x 10 4 events by year. At Tevatron, assuming 
the the running time 10 7 s and design luminosity Arevatron = 2 x 10 32 cm _2 s , we have for charm 2 — 6 x 10 9 and 
for bottom 2 — 6 x 10 7 events/year. The LHC produces the greatest rates (£lhc = 10 34 cm _2 s _1 ), giving for charm 
3.5 — 17 x 10 10 and for bottom 5.5 — 24 x 10 8 events/year. Notice the large rate for bottom at LHC. 

Finally, lets discuss the experimental separation of this reaction channel. As emphasized in Ref. [8|, although the 
photoproduction cross section to be a small fraction of the hadronic cross section, the separation of this channel is 
feasible if we impose the presence of a rapidity gap in the final state. It occurs due to the proton which is the photon 
emitter remains intact in the process. Similarly to the J/^S case, we expect that a cut in the transverse momentum of 
the pair could eliminate most part of the contribution associated to the hadroproduction of heavy quarks. Moreover, 
in comparison with the hadroproduction of heavy quarks, the event multiplicity for photoproduction interactions is 
lower, which implies that it may be used as a separation factor between these processes. As stated in Ref. pL3| . 
one way to select photoproduction events is to eliminate events where both protons breakup. This should eliminate 
almost all of the hadroproduction events while retaining most of the photoproduction interactions. In Ref. fl3j ] the 
rejection factor R, which is the probability of finding a rapidity gap with width y in a pp collision has been calculated. 
For photoproduction, the authors have obtained R = 0.04(0.005) at RHIC (LHC), requiring a single rapidity gap 
and y = 2 (For details see Section VI from Ref. ^i|). These estimates can be directly applied in our analysis. As 
photoproduction always leads to a rapidity gap, this requiriment should reject relatively few signal events, leading 
to a good signal to noise ratio for selecting these events. An important background which it should be analyzed is 
the diffractive heavy quark production in the single diffraction process |46| . We postpone this study for a future 
publication. 

In summary, we have computed the cross sections for photoproduction of open heavy quarks in pp and pp collisions. 
This has been performed using well established QCD approaches, namely the collinear and semihard factorization 
formalisms as well as saturation models within the color dipole approach. For the first time, quantitative predictions 
for these approaches are presented. The obtained values are shown to be sizeable at the current accelerators energies 
(RHIC and Tevatron) and are increasingly larger at LHC. The feasibility of detection of these reactions is encouraging, 
since their experimental signature should be suitably clear. Furthermore, they enable to constraint already in the 
current colliders the QCD dynamics since the main features from photon-proton collisions hold in proton-proton 
collisions. 
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QQ 


SAT-MOD 


SEMIHARD 


COLLINEAR 


CGC 


RHIC 


cc 


377 nb 


687 nb 


782 nb 


492 nb 




bb 


3.7 nb 


3.6 nb 


4.3 nb 


4.2 nb 


Tevatron 


cc 


1.04 fib 


2.77 fib 


3.21 fib 


1.36 fib 




bb 


13.3 nb 


24.0 nb 


29.2 nb 


19.8 nb 


LHC 


cc 


3.54 fib 


14.2 fib 


16.7 fib 


4.37 fib 




bb 


55.0 nb 


182 nb 


236 nb 


108 nb 



TABLE I: The integrated cross section for the photo-production of heavy quarks in pp(p) collisions at RHIC, Tevatron and LHC. 




W[GeV] W[GeV] 



FIG. 1: The total photoproduction cross section for charm (left panel) and bottom (right panel) . The experimental measurements 
are from DESY -HERA. 
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FIG. 3: The rapidity distribution for open charm and bottom photoproduction on pp reactions at RHIC energy \/Snn = 
500 GeV. Different curves correspond to distinct high energy QCD approaches (see text). 
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FIG. 4: The rapidity distribution for open charm and bottom photoproduction on pp reactions at Tevatron energy \J Snn = 
1.96 TeV. Different curves correspond to distinct high energy QCD approaches (see text). 




FIG. 5: The rapidity distribution for open charm and bottom photoproduction on pp reactions at LHC (CMS and/or ATLAS) 
energy t/Snn = 14 TeV. Different curves correspond to distinct high energy approaches (see text). 



